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X 



It is found that in magnetized electrostatic plasma flows the 
velocity shear couples ion-acoustic waves with ion-cyclotron 
waves and leads, under favorable conditions, to their effi- 
cient reciprocal transformations. It is shown that in a two- 
dimensional setup this coupling has a remarkable feature: it is 
governed by equations that are exactly similar to the ones de- 
scribing coupling of sound waves with internal gravity waves 
[Rogava & Mahajan: Phys. Rev. E 55, 1185 (1997)] in neu- 
tral fluid flows. Using another noteworthy quantum mechan- 
ical analogy we calculate transformation coefficients and give 
fully analytic, quantitative description of the coupling effi- 
ciency for flows with low shearing rates. 

52.30.-q, 52.35.-g, 52.35.Fp, 94.20.Bb 



I. INTRODUCTION 

It is well-known that in flows with spatially inhomo- 
geneous velocities, due to the non-self- adjointness ("non- 
normality" [1]) of their linear dynamics, collective phe- 
nomena are strongly influenced by the presence of the 
velocity shear [2]. In non-uniformly flowing plasmas, 
for example, characteristic nonperiodic modes — "Lang- 
muir vortices" [3], "ion-acoustic vortices" [4] and "dust- 
acoustic vortices" [5] — are arising due to the presence 
of the inhomogeneous velocity field. They are analogous 
to so called "Kelvin modes" [6] (originally disclosed by 
Lord Kelvin [7]) in plane Couette flows. Moreover, "tra- 
ditional" modes, present in plasmas also in the absence 
of the flow, are also interacting with "ambient" flows and 
strongly modified by the velocity shear. The modes may 
undergo transient (algebraic) growth [6], and/or linear 
(adiabatic) amplification [8] of their energy, asymptoti- 
cally persistent shear instabilities may also arise [9] due 
to the kinematic complexity of "parent" flows. 

Another important and interesting issue, related with 
the propagation of waves in flows, is related with the 
ability of the velocity shear to couple waves [10] and 
to ensure their mutual transformations. Shear-induced 
wave couplings exist in hydrodynamic systems (coupling 
of sound waves and internal gravity waves [11]), in MHD 
(both electron-proton and electron-positron plasmas), 
where velocity shear couples all three MHD wave modes 



(Alfven waves, slow magnetosonic waves and fast magne- 
tosonic waves) [12], and in dusty plasmas [5]. This phe- 
nomenon is important in a wide variety of astrophysical 
applications, including solar atmosphere and solar wind 
[13], pulsar magnetosphere [14], galactic gaseous disks 
[15], self-gravitating, differentially rotating dusty plasma 
clouds [16] and astrophysical jets [17]. 

The mission of this paper is threefold. First, it demon- 
strates that the shear-induced transformations in plas- 
mas are not characteristic only to MHD systems but they 
exist and are efficient also for low-frequency electrostatic 
waves: the ion-sound waves (ISW) and the ion-cyclotron 
waves (ICW). Second, using recently developed method 
of the evaluation of the transformation coefficients [18], 
we give full quantitative and analytic description of the 
coupling efficiency for the case of low shearing rate flows. 
Third, we point out that the mathematical description of 
this particular sort of wave coupling is exactly analogous 
to the one in neutral fluid flows: coupling of sound waves 
with internal gravity waves [11]! This remarkable cir- 
cumstance once more highlights the universality of shear- 
induced phenomena in the physics of continuous media. 



II. MAIN CONSIDERATION 

Let us consider the simplest possible setup: the equi- 
librium magnetic field Bq is homogeneous and directed 
along the X axis Bo = (So, 0,0); the mean velocity 
Uo = (Ay, 0, 0) is also directed along the X axis and has 
a linear shear along the Y axis. The equilibrium densi- 
ties of both electrons M e and ions Mi are homogeneous. 
The temperature of electrons T e is constant, while ions 
are supposed to be cold Tj = 0. 

It is well known that ISW and ICW with phase ve- 
locities considerably less than Alfven speed uj/Ux-^Va = 
Bo/ '\JAirMMi {M is the ion mass) do not perturb the 
background magnetic field, and a good approximation 
for the electric field [19] is: E = —Vip. 

Another common assumption for low-frequency elec- 
trostatic waves is that the electrons "thermalize along 
the field lines" and the electron number density can be 
described by the Boltzmann distribution [19]: 



M e = Moexp[eip/T e ]faNo[l + etp/T e 



(1) 



1 



while dynamics of ions is governed by fluid equations of 
continuity and motion: 



d t N i + V-(N i V i )=0, 



(2) 



[dt + (Vi-V)]Vi 



M 



>+-[V t xB]|. (3) 



The basic system of linearized equations for ions, de- 
scribing the evolution of small-scale, 3D perturbations in 
this flow, takes the form: 



Vtiii + Af [d x u x + dyUy + d z u z ] = 0, 



V t u x + Au y = -(e/M)d x tp, 



V t u y = -(e/M)d y ip + tt c u z , 



T> t u z 



-(e/M)d z (p - tt c ii 



v 



(4) 



(5a) 



(56) 



(5c) 



while the Poisson equation relates the perturbation of 
the electric potential ip with the number density pertur- 
bations for electrons n e and ions n>i, respectively: 



[d 2 x + d y ] + d 2 ]tp = 47re(n e -n l 



(6) 



Here fi c = eBo/Mc is the ion-cyclotron frequency and 
T> t = dt + Ayd x is the convective derivative operator. 

Employing for all perturbational variables appear- 
ing in the above equations the ansatz F(x, y, z; t) — 
F{t)exp[i{k x x + k y (t)y + k z z)\ with k y (t) = k y (0) — Atk x . 
This ansatz guarantees that V t F = exp[i(k x x + k y {t)y + 
k z z)]d t F and we can reduce the initial set of partial 
differential equations to the set of ordinary differential 
equations for the amplitudes F (see for details, e.g., [1]). 
It is convenient to write these equations in dimension- 
less notation. We define uj c = il c /C s k x as the normal- 
ized ion-cyclotron frequency and £ = Xokx, using for 
normalization conventional definitions of the ion-sound 
speed C s = (T e /M) x / 2 and electron Debye length = 
(r e /47r7v r e 2 ) 1/2 . Other notation arc: R = A/C s k x , 
t = C s k x t, p = k y /k x , (3{t) = fa - Rt, 70 = k z /k Xl 
v XiVtZ = u x , y /C s , D = i(ni/Af ), $ = ieip/MC 2 . In the 
language of these terms the set of first order, ordinary 
differential equations, derived from (4) and (5) takes the 
following form: 



D {1) = v x + p{r)vy + 1Vz , 



(7) 



Additionally we have the algebraic relation between D 
and $, which follows from the Poisson equation (6): 



D=[l+e(l+P 2 (r)+l 2 )]F 1 



(9) 



Note that (_lu c = fl c /Q p = V A jc [f7 p =(47rA^ e 2 /M) 1 / 2 
is the ion plasma frequency]. When the Debye length is 
sufficiently small (£<Cl) the oscillations can be considered 
quasi-neutral (DmF). In the forthcoming consideration 
we shall employ this approximation. However, we should 
bear in mind that due to the "k(t)-drift" , induced by 
the existence of the non-uniform motion, /3 2 (t) sooner 
or later becomes large enough and violates the quasineu- 
trality condition. 

As regards oj c , its value may be as less as greater 
than unity, depending on the relative smallness of £ and 
Va/c [Q c /flp] parameters. 

The simplest sort of waves, which may be considered in 
this system [19] are ones propagating in the XOY plane 
(7 = 0). In this case if we introduce an auxiliary nota- 
tion Y= — u> c D — P(t)v z , the system (7-8) reduces to the 
following pair of coupled second-order differential equa- 
tions: 

Y^ + ^Y + C 12 (t)v z =0, (10a) 



+ lu 2 (t)v z + C 21 (t)Y = 0, (106) 

describing coupled oscillations with two degrees of free- 
dom, with: C\ 2 {t) — C2i(r) = /3(t) as the coupling coef- 
ficients and with uoi = 1 and uj 2 (t) = uj 2 + P 2 (t) as the 
two respective eigenfrequencies. The presence of shear in 
the flow (R=/=0) makes coefficients variable and opens the 
door for mutual ICW-ISW transformations. 

In the "shcarlcss" (R = 0) limit (10) describe two in- 
dependent oscillations with frequencies: 



n 2 - 

Si 1 2 — 



1 



cc 2 +c 2 ± v /(c 2 -a, 2 2 )2 + 4C 1 2 2 



u 2 z + 1 + P 2 ± V(/3 2 +w c 2 + l) 2 -4o;2 



(11) 



that can be readily identified as ISW and ICW frequen- 
cies respectively. Corresponding eigenf unctions (some- 
times called normal variables) are: 



*i = 



(n 2 -oj 2 )Y + C 12 v z 



(12a) 



= - Rv n 



v y 1] = -P{t)<$> + u c v Zi 



(8a) 



(86) 



(fl 2 - uj 2 )v z - C 12 Y 



(126) 



,(1) 



-7$ - U> C Vy. 



(8c) 
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A. Transformation coefficient 

The coupled oscillator systems similar to (10) with 
slowly varying coefficients are well known in different 
branches of physics. Mathematical methods for their 
analysis were first developed for quantum mechanical 
problems: non-clastic atomic collisions [20] and non- 
adiabatic transitions in two level quantum systems [21]. 
Later, similar asymptotic methods were successfully ap- 
plied to various other problems [22]. 

In [18] these efficient mathematical tools were used, 
for the first time, for the study of the velocity shear in- 
duced coupling and transformation of MHD waves. The 
problem which we arc studying now, also allows thor- 
ough probing by means of this asymptotic method. We 
consider (most interesting for practical applications) case 
R <C 1, when coefficient in (10) are slowly varying func- 
tions of t and, therefore, WKB approximation is valid 
everywhere except nearby the turning points (f2j(rj) = 0) 
and resonant points [fli(r r ) = £l 2 (T r )). Using (11) one 
can check that the condition 



a (1) « fi? 



(13) 



is satisfied for both wave modes at any moment of time, 
or equivalently, none of the turning points are located 
near the real r-axis. From physical point of view this 
means, that there are no (over)reflcction phenomena [18] 
and only the resonant coupling between different waves 
modes with the same sign of 
phase velocity can occur. 

From (11) we also learn that there are two pairs of 
complex conjugated resonant points of the first order 1 : 

/3(r 1 ± ) = ± i (^ c -l), 0(-r^) = ±i(w c + l). (14) 

Therefore all the resonant points are located on the axis 
Re[/3(r)] = in the complex r-plane and consequently, 
the resonant coupling can take place only in a vicinity of 
the point where /3(r*) = 0. Generally, the time scale 
of resonant coupling At is of the order of At <~ Rr 1 / 2 
[18] and the evolution of the waves is adiabatic when 



» R 1 ' 2 . 



(15) 



If this condition is met the temporal evolution of the 
waves is described by the standard WKB solutions: 



±i j n»(r)dr 



(16) 



where D i are WKB amplitudes of the wave modes with 
positive and negative phase velocity along the x-axis, re- 
spectively. All the physical quantities can be easily found 



1 The resonant point is said to be of the order n if (fi? — fi|) 
(r — r r ) n / 2 in the neighborhood of the r r . 



by combining (12). One can check that the energies of 
the involved wave modes satisfy the standard adiabatic 
evolution condition: 



E i = Sl i (T)(\D+\' + \D. 



(17) 



Let us assume that initially [3(0) » R}l 2 and, there- 
fore, evolution of the waves is originally adiabatic. Due 
to the linear drift in the k-space, (3(t) decreases and 
when the condition (15) fails to work, the mode dynam- 
ics becomes non-adiabatic due to the resonant coupling 
between the modes. The duration of the non-adiabatic 
evolution is given by At <~ Rr x l 2 . Afterwards, when 
(3(t) <C — -R -1 / 2 , the evolution becomes adiabatic again. 
Denoting WKB amplitudes of the wave modes before and 
after the coupling region (i.e., for the t < t* — At/2 
and t > t* + At/2) by Df B and Df A respectively and 
employing the formal analogy with the S-matrix of the 
scattering theory [23] and the transition matrix from the 
theory of multi- level quantum systems [24] , one can con- 
nect Df A with Df B via the 4x4 transition matrix: 



D 



A _ 



T T± 

r|-i=p rj-i# 



(18) 



where and D^; while 1x2 matrices and T, its Her- 
mitian conjugated matrix T*, T ± , and T T are 2x2 
matrices. 

None of the turning points are located near the real r- 
axis and, therefore, only wave modes with the same sign 
of the phase velocity along the x-axis can effectively in- 
teract. It is well known [24] that in this case components 
of and T T are exponentially small with respect to 
the large parameter l/R and can be neglected. Conse- 
quently, (18) decomposes and reduces to: 



TD i 



T*D^ 



(19a) 



(196) 



Since all coefficients in the (12) are real and C\2 = C21, 
the matrix T is unitary [25], and 



Ei 



= 1. 



(20) 



Generally, this equation represents conservation of the 
wave action. When R -C 1 it transcribes into the en- 
ergy conservation throughout the resonant coupling of 
wave modes [18]. The components of the transition ma- 
trix in (18) are complex, i.e., the coupling of different 
wave modes changes not only the absolute values of Df , 
but also their phases. The value of the quantity |Ti 2 | 2 
represents a part of the energy transformed during the 
resonant coupling of the modes. The absolute values 
of the transition matrix components |Ti2 1 and IT21I are 
called the transformation coefficients of corresponding 
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wave modes. Unitarity of the T provides an important 
symmetry property of the transition matrix: 



e (2) + [w 2 + 1 (r)] e + /3(r)V - 0. 



(266) 



\T, 



(21) 



i.e., transformation coefficients are reciprocally equal to 
each other. 

It is well known, that if in the neighborhood of the 
real r-axis only a pair of complex conjugated first or- 
der resonant points r r and r* exists, the transformation 
coefficient is [18]: 



-12 



cxp 



Im / (fii - n 2 )dr 



[1 + OiR 1 ' 2 )]. 



(22) 

As it was mentioned earlier, coupling between ICW 
and ISW can be effective if \u c — 1| <C 1. Otherwise 
transformation coefficient is exponentially small with re- 
spect to the large parameter R^ 1 [24,25]. If this condition 
is satisfied the resonant points [see (14)] tend to the 
real r-axis and, hence, the effective coupling is possible. 
Noting, that according to (13), in the neighborhood of 
the resonant points: 



Ox - 2 w V/3 2 (t) + (w c -1) 2 



(23) 



and then from (22) one can readily obtain: 

|T 12 |«exp[-^ c -l)74.R]. (24) 

From the (14) we also see that if uj c —* 1, the resonant 
points tend to the real r-axis. Then from (24) it follows 
that |Ti2 1 — * 1, i.e., one wave mode is totally transformed 
into the another. On the other hand, if \u) c — 1| 3> R 1 ^ 2 , 
transformation is negligible. 

In the case of moderate or high shearing rates, similarly 
to the case R <C 1, one can show that the WKB approx- 
imation is valid only when /3(r) 3> 1,R. It means that 
the asymptotic problem can still be formulated. How- 
ever, when the shearing rate is not small, non-adiabatic 
evolution of the modes consists of both transformation 
and reflection phenomena. From mathematical point of 
view it means that, in general, all the components of 
the transition matrix significantly differ from zero. Con- 
servation of the wave action remains valid and provides 
following important relation [18]: 



(25) 



However, it should be also stressed that for high shear- 
ing rates no analytical expressions for the components of 
the transition matrix can be written explicitly. 



B. Hydrodynamical analogy 

The remarkable feature of our governing equations (10) 
is that they are almost identical with the pair of equations 
from [11] (numbered there as equations (16) and (17)): 



V> (2) + V + 0(r)e = 0, 



(26a) 



These equations also describe coupling of two wave 
modes. But this is totally different physical system: 
shear flow of a gravitationally stratified neutral fluid, 
which sustains sound waves and internal gravity waves. 
In [11] it was shown that these modes are coupled through 
the agency of the shear and may effectively transform into 
each other, providing the condition VF~1 is met. 

The presence of this analogy implies that all the de- 
tails of the transformation coefficient asymptotic analy- 
sis, which were given above, can also be applied to the 
named hydrodynamic example of the shear induced wave 
couplings. The only factual difference is that we have to 
replace the dimensionless frequency of electrostatic ion- 
cyclotron waves llJ c by the characteristic dimensionless 
frequency W of the internal gravity waves. 

This analogy allows also to predict that within the elec- 
trostatic problem we should have a new kind of electro- 
static beat waves just as they exist in the hydrodynamical 
problem! Beat waves are excited when an initial pertur- 
bation propagates almost along the flow axis (/3o<d and 
when, additionally, r* is sufficiently large in comparison 
with the beat period flb=fli — £1 2 - 



III. CONCLUSION 

Summarizing main properties of the resonant coupling 
of ICW and ISW for small shearing rates (R <C 1) we can 
state that: (a) Only the wave modes with the same sign 
of the phase velocity can effectively interact - there are 
no reflection phenomena; (b)The duration of the effective 
coupling of the modes is of the order of Ar ~ R~ x / 2 , i.e., 
resonant coupling is slow compared to the wave period 
m <~ 1 but fast enough compared to the adiabatic change 
of the system parameters r a ~ (c) The total energy 
of the modes is conserved during the resonant coupling - 
the transformed wave is generated at the expense of the 
energy of the initial wave mode; (d) The mode coupling 
process is symmetric - transformation coefficient of one 
mode into another one equals the coefficient of the inverse 
process, (e) The transformation coefficients are given by 
(22) and (24). 

None of these features remain valid for moderate and 
high shearing rates. If R <C 1 is not satisfied, there are 
no 'long' and 'shirt' timescales in the problem and all 
the processes have approximately the same characteristic 
timescale. Hence, the coupling of the modes represents 
some mixture of transformation and reflection precesses, 
that are accompanied by the energy exchange between 
the waves and the background flow. 

The discovered ICW-ISW transformations are likely 
to be important in a number of applications. One possi- 
ble example is the problem of the ICW observed by low 
altitude satellites and ground based magnetometers [26] . 
Observational surveys indicate that these waves are cor- 
related with the ICW observed in the equatorial magne- 
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tosphcrc. However, theoretical studies of the ICW prop- 
agation from the magnetosphere to the ground suggest 
that these waves can not penetrate through the Buchs- 
baum resonance and can not reach ionospheric layers of 
the atmosphere. Thus, one could expect that the mag- 
nctosphcric ICW shouldn't be correlated with the iono- 
spheric ICW, while observational evidence shows the cor- 
relation. Recently Johnson and Cheng [26] reconsidered 
this problem and found that strong mode coupling oc- 
curs near the He + and + resonance locations. They 
argued that this coupling may help the equatorial ICW 
to penetrate to ionospheric altitudes. 

It seems reasonable to admit that the velocity- shear- 
induced ICW-ISW coupling may provide yet another 
mode transformation mechanism, which in conjunction 
with the one, found by Johnson and Cheng, may account 
for the penetration of the ICW through the Buchsbaum 
resonance. This may work in a quite similar way to the 
scenario given in the [27] for the penetration of the fast 
magnetosonic waves (FMW) from the chromosphere to 
the corona through the solar transition region. The idea 
is that a fraction of the FMW transforms into the Alfven 
waves (AW). The latter go through the transition region 
up to the solar corona, where they again become trans- 
formed into the FMW. In this way 'shear-induced oscilla- 
tions of solar MHD waves' [27], ensures substantial trans- 
port of the FMW through the transition region. It seems 
plausible to admit that the similar reciprocal 'swinging' 
of ICW and ISW may allow some fraction of the ICW 
to penetrate through the Buchsbaum resonance and to 
reach the low ionospheric altitudes. 

Finally, the remarkable exact analogy of the ICW-ISW 
coupling with the coupling of internal gravity waves and 
sound waves in hydrodynamic flows [11] points out, once 
again, at the universal character of the velocity shear 
induced phenomena in the physics of fluids and plasmas. 

IV. ACKNOWLEDGEMENTS 

Andria Rogava wishes to thank Abdus Salam Interna- 
tional Centre for Theoretical Physics for supporting him, 
in part, through a Senior Associate Membership Award. 



[1] L. N. Trefethen, A. E. Trefethen, S. C. Reddy, and T. A. 

Driscoll, Science 261, 578 (1993); W.O. Criminale and P.G. 

Drazin Stud. Appl. Maths., 83, 123 (1990). 
[2] F. Volponi, Z. Yoshida, and T. Tatsuno, Phys. Plasmas 7, 

4863 (2000). 

[3] A. D. Rogava, G. D. Chagelishvili, and S. M. Mahajan, 
Phys. Rev. 55, 7103 (1998). 

[4] A. D. Rogava, G. D. Chagelishvili, and V. I. Berezhiani, 
Phys. Plasmas 4, 4201 (1997); 



[5] S. Poedts, G. R. Khujadze, and A. D. Rogava, Phys. Plas- 
mas 7, 3204 (2000) 

[6] P. Marcus and W. H. Press, J. Fluid Mech., 79, 525 (1977). 
[7] Lord Kelvin (W. Thomson), Phil. Mag. 24, Ser. 5, 188 
(1887); 

[8] M. Butler and B. F. Farrell, Phys. Fluids (A) 4, 1637 
(1992); G. D. Chagelishvili, A. D. Rogava, and I. N. Segal, 
Phys. Rev. (E) 50, R4283 (1994); G. D. Chagelishvili, G. 
R. Khujadze, J. G. Lominadze, and A. D. Rogava, Phys. 
Fluids 9, 1955 (1997). 

[9] S. M. Mahajan, and A. D. Rogava, Ap.J., 518, 814 (1999). 
[10] G. D. Chagelishvili, A. D. Rogava, and D. G. Tsiklauri, 
Phys. Rev. (E) 53, 6028 (1996). 

[11] A. D. Rogava and S. M. Mahajan, Phys. Rev. (E) 55, 
1185 (1997). 

[12] A. D. Rogava, S. M. Mahajan, and V. I. Berezhiani, Phys. 
Plasmas 3, 3545 (1996); G. D. Chagelishvili, A. D. Ro- 
gava, and D. G. Tsiklauri, Phys. Plasmas 4, 1182 (1997); T. 
Tatsuno, F. Volponi, and Z. Yoshida, Phys. Plasmas 8, 399 
(2001). 

[13] S. Poedts, A. D. Rogava, and S. M. Mahajan, Ap. J. 505, 
369 (1998); A.D. Rogava, S. Poedts, and S. M. Mahajan, 
Astron. Astrophys. 354, 749 (2000). 

[14] S. M. Mahajan, G. Z. Machabeli, and A. D. Rogava, Ap. 
J. Lett 479, 129 (1997). 

[15] Z. H. Fan, and Y. Q. Lou, Mon. Not. Royal. Astron. Soc, 
291, 91 (1997); A. D. Rogava, S. Poedts, and S. Heirman, 
Mon. Not. Royal Astron. Soc, 307, L31 (1999); S. Poedts 
and A. D. Rogava, Astron. Astrophys. 385, 32 (2002). 

[16] A. D. Rogava, S. Poedts, and Z. N. Osmanov, Phys. Plas- 
mas 11, 1655 (2004). 

[17] A.D. Rogava, S. M. Mahajan, G. Bodo, and S. Massaglia, 
Astron. Astrophys. 399, 421 (2003); A. D. Rogava, G. Bodo, 
S. Massaglia, and Z. Osmanov, Astron. Astrophys. 408, 401 
(2003). 

[18] G. Gogoberidze, G. D. Chagelishvili, R. Z. Sagdeev and 
J. G. Lominadze, Phys. Plasmas, 11, 4672, (2004). 

[19] F. F. Chen, Introduction to Plasma Physics (Plenum 
Press, New York and London, 1977). 

[20] E. C. J. Stuekelberg, Helv. Phys. Acta. 5, 369 (1932). 

[21] C. Zener, Proc. R. Soc. London Ser. A 137, 696 (1932); 
L. D. Landau, Phys. Z. Sowjetunion 2, 46 (1932). 

[22] D. G. Swanson, Theory of Mode Conversion and Tun- 
nelling in Inhomogeneous Plasmas (John Wiley and Sons, 
New York, 1998). 

[23] T. Kopaleishvili, Collision theory: a short course (World 
Scientific Publishing Corporation, 1995), p. 38. 

[24] L. D. Landau and E. M. Lifschitz, Quantum Mechanics 
(Non-Relativistic Theory) (Pergamon Press, Oxford, Eng- 
land, 1977), p. 304. 

[25] M. V. Fedoriuk, Asymptotic methods for ordinary differ- 
ential equations (Nauka, Moscow, 1983), p. 324. 

[26] J. R. Johnson, and C. Z. Cheng, Geophys. Res. Lett., 26, 
671 (1999). 

[27] S. Poedts, A. D. Rogava, and S. M. Mahajan, Space Sci. 
Rev., 87, 295 (1999). 



•5 



